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ABSTRACT 

Flaviviral RNA-dependent RNA polymerases 
(RdRps) initiate replication of the single-stranded 
RNA genome in the absence of a primer. The 
template sequence 5'-CU-3' at the 3'-end of the 
flaviviral genome is highly conserved. Surprisingly, 
flaviviral RdRps require high concentrations of the 
second incoming nucleotide GTP to catalyze de 
novo template-dependent RNA synthesis. We show 
that GTP stimulates de novo RNA synthesis by RdRp 
from Japanese encephalitis virus (jRdRp) also. 
Crystal structures of jRdRp complexed with GTP 
and ATP provide a basis for specific recognition of 
GTP. Comparison of the jRdRpcTP structure with 
other viral RdRp-GTP structures shows that GTP 
binds jRdRp in a novel conformation. Apo-jRdRp 
structure suggests that the conserved motif F of 
jRdRp occupies multiple conformations in absence 
of GTP. Motif F becomes ordered on GTP binding 
and occludes the nucleotide triphosphate entry 
tunnel. Mutational analysis of key residues that 
interact with GTP evinces that the jRdRpcTP struc- 
ture represents a novel pre-initiation state. Also, 
binding studies show that GTP binding reduces 
affinity of RdRp for RNA, but the presence of the 
catalytic Mn^^ ion abolishes this inhibition. 
Collectively, these observations suggest that the 
observed pre-initiation state may serve as a check- 
point to prevent erroneous template-independent 
RNA synthesis by jRdRp during initiation. 

INTRODUCTION 

Japanese encephalitis virus (JEV) is a positive-strand 
RNA virus that belongs to the genus Flavivirus of the 



family Flaviviridae. Other members of this genus include 
dengue (DENY), West Nile (WNV), yellow fever and tick- 
borne encephahtis viruses (1). JEV is a neurotropic virus 
and the leading cause of viral encephalitis. The disease 
caused by this virus is endemic in China, South and 
South-east Asia. Conservative estimates state that 
~50000 cases are reported annually, one-third of which 
are fatal and around half result in permanent neurological 
damage in survivors (2). JEV is thought to be responsible 
for hundreds of fatalities — mostly children — that occur 
every year in northern India due to viral encephahtis. 
Although vaccines are available, a rigorously tested effect- 
ive low-cost vaccine that will afford long-term protection 
against infection is still lacking. As in the case of other 
pathogenic tlaviviruses, therapeutic agents to combat JEV 
infection and to prevent disease progression are not 
available. 

The JEV genome is represented by a single-stranded 
RNA (ssRNA) molecule ~ll-kb long, with a 5' cap 
(m7G5'ppp5'A) and no polyadenylated tail. This viral 
genome codes for three structural (C, prM and E) and 
seven non-structural (NSl, NS2A, NS2B, NS3, NS4A, 
NS4B and NS5) proteins (3). The primary translation 
product of the viral genome is a 370-kDa polyprotein, 
which is processed by a host signalase and viral protease 
to yield structural and non-structural proteins. 

NS5 is the largest (~105kDa) and the most conserved 
protein among the members of the Flavivirus genus (1,4). 
The N-terminal region of this protein (~30 kDa) houses a 
methyltransferase activity, which is responsible for adding 
the methyl group to the 5' RNA cap structure (5,6). 
Toward the C-terminus is present the RNA-dependent 
RNA polymerase (RdRp) module (~70kDa) that is 
critical for replication of the viral genome (7,8). NS5 as- 
sociates strongly with NS3 and along with other viral 
proteins, host factors and genomic RNA "forms" a repli- 
cation complex (RC) (9,10). Initially, the negative strand is 
synthesized to give rise to a double-stranded replicative 
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form (dsRF) within the membrane-associated RC. The 
RC is localized inside double-layered membrane vesicles, 
presumably to protect the dsRF from host anti-viral 
factors (11,12). The RF is then used as a template to 
generate replicative intermediates. These replicative inter- 
mediates ultimately resolve to give rise to single-stranded 
genomic RNA and dsRF on completion of one round of 
synthesis. Thus, the synthesis of progeny RNA occurs in a 
semi-conservative and asymmetric manner on recycling 
dsRF (13). In addition to its core function in viral 
genome replication, NS5 has been implicated in regulating 
expression of host genes involved in anti-viral responses 
(4,14). 

Unlike replicative DNA polymerases, viral RdRps use 
diverse strategies to initiate RNA synthesis (15). The poly- 
merase of Picornaviridae family uses a protein primer (16), 
while influenza virus of the Orthomyxoviridae family can 
use an unusual cap snatching method to prime rephcation 
(17). For members of the families Flaviviridae, Cystoviridae 
(phi6 bacteriophage) and Reoviridae (rotaviruses), the 
RdRp enzyme can initiate RNA synthesis at the template 
termini de novo without the need for any primer (15). Once a 
phosphodiester bond is formed between the first two nucleo- 
tides, the enzyme switches from initiation to elongation 
mode and mediates template-dependent extension of the 
dinucleotide primer (18). 

In the case of DENV RdRp, a high concentration of 
GTP is required for initiation of synthesis of RNA (18- 
20). The presence of the CU sequence at the 3'-end of the 
genomic RNA imphes that GTP should be the second 
incoming nucleotide. The CU sequence is conserved in 
the genomes of most flaviviruses including JEV, at 
3'-ends of both positive and negative strands. The require- 
ment for higher concentration of GTP is also true for 
RdRp from distant viruses like hepatitis C virus (HCV; 
genus Hepacivirus; family Flaviviridae) (21), classical 
swine fever virus (genus Pestivirus; family Flaviviridae) 
(22) and phi6 bacteriophage (family Cystoviridae) (23). 
It should be mentioned here that the genomes of these 
viruses exhibit distinct terminal dinucleotide sequences 
that can be different for the negative and positive strands. 

We show here that JEV RdRp (jRdRp) requires high 
concentration of GTP for initiation, but not for elong- 
ation. We have determined two crystal structures of 
jRdRp in complex with GTP (2.4 A) and ATP (2.3 A). 
In addition, we have obtained a structure of RdRp in 
the apo-state (3.65 A). The structures of the jRdRpoxp 
and jRdRpATP complexes provide a basis for specific rec- 
ognition of GTP by this enzyme. Additionally, these struc- 
tures showed that the conserved motif F undergoes 
substantial ordering and occludes the nucleotide triphos- 
phate (NTP) entry tunnel on GTP binding. Comparison 
of the jRdRpoTP structure with that of other viral RdRp- 
NTP complexes showed that the GTP molecule is bound 
in a novel conformation. A number of point mutants of 
jRdRp were prepared on the basis of the comparison of 
the jRdRpoxp structure with that of full initiation 
complex of RdRp from the bacteriophage phi6. These 
mutants were assessed for their ability to catalyze initi- 
ation and elongation of RNA synthesis. These experi- 
ments showed that jRdRpoTP structure represents a 



novel pre-initiation state of genome replication. 
Additionally, it was seen that GTP inhibits binding of 
jRdRp to RNA and that the presence of Mn^^ ions 
abolishes this inhibitory effect. Also, this inhibitory 
effect of GTP on RNA binding is lost or reduced in a 
mutant of RdRp that shows reduced initiation activity. 
Taken together, these studies allow formulation of a 
possible sequence of events during initiation that 
prevents non-templated and wasteful RNA synthesis 
that could be erroneous. 



MATERIALS AND METHODS 

Cloning, expression and purification of the jRdRp enzyme 

The gene segment corresponding to the RdRp module 
(corresponding to amino acid 272-905 of the NS5 
protein) was amplified from pGEX6Pl plasmid containing 
the NS5 gene of JEV strain P20778 (GenBank accession 
number AF080251). The amplified product was purified 
and hgated into the pGEX6Pl that would yield a fusion 
polypeptide with a GST tag at the N-terminus connected 
to the RdRp module by a hnker with the PreScission 
protease site. Sequencing was carried out for a few 
clones and one without any mutations was selected for 
expression and purification of the RdRp protein. For 
making site-specific mutations in the jRdRp enzyme, 
QuikChange Lightning Multi Site-Directed Mutagenesis 
Kit (Agilent Technologies) was used. Primers were 
designed, and mutants were prepared according to manu- 
facturer's instructions. 

The GST-jRdRp fusion protein was expressed in 
Escherichia coli C41(DE3) cells. Five liters of LB 
medium (100|ig/ml ampicilhn) was inoculated with 25 ml 
of starter culture and incubated at 37°C at 200 rpm. When 
the optical density at 600 nm reached 0.6-0.8, isopropyl 1- 
thio-P-D-galactopyranoside was added to the culture at a 
final concentration of 0.25 mM and was incubated for 18 h 
at 18°C. Cells were harvested by centrifugation, resus- 
pended in 25 mM Tris-Cl buffer, pH 8.0 (4°C), 500mM 
NaCl, 5% glycerol, 2mM dithiothreitol (DTT), 0.01% 
IGEPAL CA-630 and 1 mM phenyl methyl sulfonyl fluor- 
ide and then frozen in — 80°C until use. Cells were lysed by 
sonication, and the lysate was clarified by centrifugation 
at 17 000 rpm for 45min at 4°C. Glutathione-sepharose 
resin pre-equilibrated with buffer A (25 mM Tris-Cl 
buffer, pH 8.0, 250 niM NaCl, 5% glycerol and 2mM 
DTT) was incubated with the clarified supernatant for 
3h. The beads were then washed with buffer B (25 mM 
Tris-Cl buffer, pH 8.0, 750mM NaCl, 5% glycerol and 
2mM DTT) to remove non-specifically bound proteins, 
and then again with buffer A. The resin was then 
incubated with PreScission protease overnight to cleave 
the GST tag. jRdRp protein was eluted with buffer A, 
concentrated and further purified by size exclusion chro- 
matography using a Superdex-200 column equilibrated 
with buffer C (25 mM Tris-Cl, pH 7.4, 250 mM NaCl, 
2% glycerol and 2mM DTT). The presence of jRdRp in 
the major peak was confirmed by sodium dodecyl 
sulphate-polyacrylamide gel electrophoresis analysis, and 
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the purified jRdRp was concentrated to 8-12mg/ml and 
stored at -80°C. 

Primer-free initiation assay 

Primer-free de novo RNA synthesis was assessed by tracking 
the incorporation of 6-carboxy Fluorescein (6FAM)-labeled 
UTP in the presence of a template. The incorporation will 
occur only if a dinucleotide primer is formed and extended 
by the RdRp enzyme. Hence, in the absence of productive 
initiation, addition of labeled UTP to the growing primer 
will not occur. This method avoids the use of radioactive 
NTPs, which constitute a significant health and environmen- 
tal hazard (24). The 13-nt RNA template (20pmole) (5'-UA 
AUUCCCCCCCU-3') was mixed with jRdRp protein 
(lOpmole) along with ATP (10 nM), CTP (10|iM), FAM- 
UTP (10 nM) and GTP (100 |iM) in a 20-|.il reaction mix. 
The buffer used for initiation was 50 mM Tris-Cl, pH 8.0 (at 
30°C), ImM DTT, 5% glycerol and 0.05mg/ml bovine 
serum albumin. To assess the role of GTP and ATP in ini- 
tiation, reactions were carried out with varying concentra- 
tions of these two NTPs with the concentration of the other 
NTPs maintained at 10|iM. The concentration of divalent 
ion Mn^^ was 2.5 mM in the reaction mix. The reaction 
mixes were incubated at 30° C for 3h, and the reactions 
were stopped by adding lOiil of a solution containing 
20niM ethylenediaminetetraacetic acid and 80% 
forniamide. The samples were heated for 5min at 98°C 
followed by cooling on ice for 5min. Next, the reaction 
products were resolved on a 20% polyacrylamide gel con- 
taining 8M urea that was prepared and run in Ix Tris- 
borate + ethylenediaminetetraacetic acid buffer. Resolved 
products on the gel were observed by excitation at 488 nm, 
and the bands were visualized and recorded using Pharos 
Imager (BioRad Corp.). The intensities of the observed 
bands were quantified using Quantity one, ID analysis 
software. The initiation activity for each reaction mixture 
was quantified by calculating the cumulative intensity of all 
the product bands in the corresponding lane on the gel. 

A similar protocol was used to assess the effect of site- 
specific mutations in the jRdRp enzyme on primer-free 
initiation. For estimating relative activity of mutants, the 
activity of wild-type (wt) was taken as 100%. 

Primer extension assay 

To probe the elongation activity of jRdRp, a primer (5'-GU 
UCACACAGAUAAACUUCU-3') with a 6-FAM labeled 
at the 5'-end was used as an alternative to radioactive 
labelling (24). RNA template (5'-AGAUCCUGCCCCCC 
CCCCCAAGAAGUUUAUCUGUGUGAAC-3') and the 
labeled primer were mixed in equimolar quantity in 1 x TE 
buffer and annealed by heating at 90°C for 5 min followed 
by cooling to room temperature (25°C). The annealed 
RNA template-primer duplex (20 picomole) was incubated 
with jRdRp (10 picomole) in 20|.d of 50 mM Tris-Cl, pH 
8.0 (37°C), 1 mM DTT, 5% glycerol and 0.05mg/ml bovine 
serum albumin. All four NTPs were added to the reaction 
volume to a final concentration of 100 |aM, and the catalytic 
ion Mn2+ (MnCh) was maintained at concentration of 
2.5 mM. The reaction mixtures were incubated at 37°C 
for 3h, and the samples were then processed using a 



protocol similar to that for initiation assay. The intensity 
of the unextended primer in each lane was quantified using 
Quantity one, ID analysis software. The values obtained 
for each reaction mixture were divided by that obtained for 
a sample wherein reaction cannot occur to calculate the 
fraction of total primer that was extended by jRdRp. To 
quantify the elongation activity of different mutants of 
jRdRp relative to wt enzyme, the fraction of total primer 
extended by wt RdRp was taken as 100%). 

Crystallization 

jRdRp (0.08 mM) was incubated with GTP/ATP (5mM) 
for 30 min before setting crystallization trials. After 
screening and optimization, the best crystals were obtained 
in 10-18% wt/vol PEG 4000, lOOmM Tris-Cl buffer, pH 
8.0 (at 25°C), and 7% v/v 1-propanol by hanging drop 
method. The size of the crystals was improved by 
macroseeding. For data collection, the crystals were 
cryoprotected by soaking for 1 min in reservoir solutions 
containing 5, 10, 15 and 20% glycerol, respectively, and 
then flash frozen in hquid nitrogen. 

Crystals of apo-jRdRp (0.096 mM) were obtained in 
10-15% PEG 5000 mono methyl ether, lOOmM Tris-Cl, 
pH 8.0 (at 25°C), and 200 mM NaSCN. The crystals were 
cryoprotected using a protocol similar to that for crystals 
of jRdRpoTP- 

Structure determination 

X-ray diffraction data were collected at the BM14 
beamhne of ESRF, and a single crystal was used for 
data collection. Data were indexed, integrated, scaled 
and merged using HKL2000 (25). 

The jRdRpoTP structure was determined by molecular 
replacement (MR) using MolRep and Phaser programs in 
CCP4i suite (26). The structure of apo- RdRp from WNV 
(2HFZ) and DENV (2J7W) was used as search models 
(27,28). For successful MR, search models corresponding 
to the pahn, fingers and thumb domains of 2HFZ were 
prepared. WNV RdRp provided the best solutions for the 
thumb and palm domains with PHASER, and the collated 
coordinates were fixed during a MOLREP run to obtain a 
solution for the fingers domain. The composite solution 
was then subjected to rigid body refinement followed by 
B-factor and positional refinement using CNS (29). 
Subsequently, electron density maps (2Fo-Fc, Fo-Fc and 
composite simulated annealing omit map) were calculated 
and showed good fit between chemical model and electron 
density. These maps showed clear electron density corres- 
ponding to the GTP molecule near the active site. GTP was 
placed in this electron density, the sequence of 2HFZ was 
slowly changed to that of jRdRp and every change in the 
model was subjected to crystallographic refinement. After 
several rounds of model building and refinement, the auto- 
matic water picking procedure of CNS Package were used 
to place water molecules. This was followed by further 
refinement, and on convergence of the i?free value, 
Translation/Libration/Screw (TLS) refinement was carried 
using Refmac program in CCP4i (30). 

The jRdRpATP structure was determined by MR with 
the jRdRpGTP structure as a search model. Before MR, 
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the coordinates corresponding to GTP and water mol- 
ecules were removed from the search model. MR with 
Phaser yielded a clear solution, and after a round of 
rigid body and b-factor refinement, electron density 
maps were calculated. These maps showed density for 
the ATP molecule in the active site, but the density for 
base moiety was broken. ATP was placed in this density, 
and this was followed by model building, refinement and 
placement of water molecules iteratively until the i^tvee 
value converged. Finally, TLS refinement was carried 
out using the Refmac program in CCP4i. 

The apo-jRdRp structure was determined by MR using 
the jRdRpoTP structure as a search model in Phaser. 
Before MR, the coordinates corresponding to GTP and 
water molecules were removed from the search model. 
Refinement was carried out in PHENIX with NCS con- 
straints until convergence of the R^^^^ value (31). Data 
collection and refinement statistics for jRdRpoxp. 
jRdRpATP and apo jRdRp are provided in Table 1. 
Coordinates and structure factors for the three structures 
have been deposited with the Protein Data Bank with the 
following accession codes: 4HDG (jRdRpoTp). 4HDH 
(jRdRpATp) and 4MTP (apo- jRdRp). 

Fluorescence anisotropy measurements to monitor binding 
of RdRp and RNA 

Fluorescence anisotropy was measured using a 
SpectraMax M5 microplate reader (Molecular Devices). 
jRdRp and ssRNA (5'- XAUAAACUUCU-3', where X 



is 6FAM label) were mixed and equilibrated for 20 min at 
room temperature in a binding buffer (50 mM Tris-Cl, pH 
8.0, at 25° C, 5% glycerol, 1 mM DTT). RNA concentra- 
tion was maintained at 5nM, and jRdRp concentration 
was varied from 0 to 4000 nM. When added, Mn^^ ion 
was present in 100-fold molar excess of the protein. To 
study the effect of GTP on jRdRp-ssRNA binding, the 
mixture of jRdRp (1 \\M) and FRIO (5 nM) was incubated 
with varying amounts of GTP. The ratio of the concen- 
tration of jRdRp to that of GTP was varied from 1:0.1 to 
1:400. Fraction of RNA bound (Fsb) was calculated using 
the following equation (32), 

P _ AoBs— Ap 

(Ab-Aobs)G+Aobs-Af' 

where Aqbs is the observed anisotropy for a particular 
ratio of jRdRp and labeled RNA. Ab and Ap denote 
anisotropics of bound and free-labeled RNA, respectively. 
Q is the ratio of fluorescence intensities of bound and free- 
labeled RNA measured in the experiment. To determine 
the apparent dissociation constants (Kj), fraction of RNA 
bound (Fsb) was plotted as a function of total RdRp con- 
centration (Rt), and the data were fit to the following HiU 
equation using non-hnear regression in Origin 8 software, 
where n is the average number of interacting sites: 

(Rt-(5 * Fsb))" 
K3+(Rt-(5*Fsb))" 



Table 1. Data collection and refinement statistics 





jRdRpGTp 


jRdRpATp 


Apo jRdRp 


Data collection ^ 








Wavelength (A) 


0.95373 


0.95373 


0.95373 


Space group 


P2i2,2 


P2,2,2 


P2,2,2, 


Unit cell dimensions-a, b, c (A) 


144.6, 86.6, 112.7 


144.8, 86.5, 112.8 


119.5, 174.0,182.7 


Angles-a, (3, j (°) 


90.00, 90.00, 90.00 


90.00, 90.00, 90.00 


90.00, 90.00, 90.00 


Resolution (A) 


50-2.38 


50-2.28 


56.78-3.65 




8.4 (47.3) 


7.0 (43.5) 


15.5 (56.9) 


I/5I 


17.0 (2.7) 


21.2 (2.8) 


6.7 (2.4) 


Completeness (%) 


99.73 (99.50) 


99.4 (97.4) 


99.6 (99.9) 


Redundancy 


5.3 (5.3) 


5.0 (4.0) 


4.8 (4.8) 


Number of reflections 


57 470 


65 396 


428 66 


Refinement 










22.0/27.6 


22.4/27.0 


23.4/27.3 


Number of atoms 








Protein 


9838 


9830 


18 702 


Nucleotide 


64 


62 




Zn"^ ion 


4 


4 


8 


Water 


455 


422 




Average B-factors (A^) 








Protein 


36.6 


38.3 


47.1 


Nucleotide triphosphate 


67.0 


64.1 




Zn^^ ion 


38.6 


36.6 


44.9 


Water 


36.6 


38.1 




R.M.S.D 








Bond lengths (A) 


0.017 


0.019 


0.002 


Bond angles (°) 


1.587 


1.567 


0.64 



Values in brackets refer to the highest resolution shell. 

''Riree was Calculated with 7% of reflections excluded from the whole refinement procedure for jRdRpcxp and jRdRpAxp- For apo jRdRp, 5% of 
reflections were excluded. 
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RESULTS 

GTP and initiation of RNA synthesis by jRdRp 

The two nucleotides CU at the 3'-ends of positive and 
negative strands are conserved in the genomes of most 
flaviviruses and, therefore, could be essential for initiation. 
Additionally, RdRp molecules are known to go through 
rounds of abortive initiation before switching to elong- 
ation mode (20). To test the role of GTP in initiation, 
we designed a template RNA (5'-UAAUUCCCC 
CCCU-3') that retains the conserved dinucleotide CU at 
the 3'-end and will not loop back to prime RNA synthesis. 
For this template, we observed maximum RNA synthesis 
at high GTP concentrations (>100|.iM) with the concen- 
tration of other nucleotides maintained at lOuM 
(Figure lA and B). At high concentrations of GTP, a 
number of products corresponding to different lengths 
of RNA were synthesized giving rise to a ladder-like 
pattern on the gel (Figure lA). The ladder-like pattern is 
probably due to polymerase slippage caused by the stretch 
of consecutive cytosines in the template, a phenomenon 
seen routinely in the case of DNA-dependent DNA 
polymerases (33). Additionally, RdRp molecules are 
known to go through rounds of abortive initiation 
before switching to elongation mode and this may 
heighten the effect of polymerase slippage (20). Unlike 
the case of GTP, there was no comparable increase in 
the amount of initiation products with increasing concen- 
trations of ATP (Figure IB and C). 

It was seen that Mn^"^ ions were necessary for initiation, 
and no products were formed in the presence of Mg^^ 
ions. The requirement for Mn^^ was also observed in 
earlier reports on biochemical characterization of recom- 
binant jRdRp (34,35). However, for replication complexes 
isolated from membranes of cells infected with JEV, it has 
been shown that in addition to low concentrations of 
Mn-+ (0.05-0.4 mM), higher concentrations of Mg^+ 
(2-10 niM) can also support RNA synthesis (36). 

We also assessed the effect of varying concentrations of 
GTP on elongation of a labeled primer annealed to a 
template. The conditions for primer extension were 
standardized, and optimal activity was obtained in the 
presence of Mn^"^ ions at 37°C. In the substrate used for 
the primer extension assay, the label is present at the 
5'-end of the primer and the ladder-hke pattern observed 
in the gel appears due to dissociation of the polymerase 
after a few rounds of replication. There was no significant 
qualitative or quantitative difference in the products 
formed in the presence of different concentrations of 
GTP or ATP (other NTPs were maintained at a concen- 
tration of 100|.iM) (Figure ID and E). Taken together, 
these experiments show that high concentration of GTP 
aids initiation but is not necessary for elongation. 

Structure of the RdRp enzyme in jRdRpcxp and 
jRdRpATP 

The two complexes jRdRpoxp and jRdRpAxp crystallized 
in the space group P2i2i2 with the cell constants of 
a = 144.6 A, b = 86.6A, c = 112.7 A and a = p = y = 90° 
(Table 1). The jRdRpoxp complex was refined to a final 



^fi-ee of 27.6% and i?cryst of 22.0% (2.4 A resolution). In 
the case of jRdRp atp, the final i^f^ee md R„yst are 27.0% 
and 22.4%, respectively (2.3 A resolution). The asymmet- 
ric unit contains two representatives of the jRdRpNxp 
complex. Each complex contains one molecule of jRdRp 
with residues from 274 to 315 and 321 to 889 (the regions 
316-320 and 890-905 are disordered) (Figure 2A), one 
molecule of GTP/ATP, water molecules and two Zn^^ 
ions. Each Zn^^ ion is present at a distinct location 
away from the active site (Supplementary Figures SI 
and S2). These locations are different from the ones 
observed for the catalytic ion (Mn^"^/Mg^"^) in structures 
of ternary complexes of other viral RdRps. Zn^"^ ions in 
jRdRp are present at the same location as in the case of 
DENY RdRp (2J7W) and WNV RdRp (2HFZ). 

The presence and conformation of the bound GTP was 
confirmed through calculation of a simulated annealing 
Fo-Fc omit map (Figure 2B). The structure of the two 
RdRp molecules in the asymmetric unit is nearly identical 
and can be superimposed with a root-mean-square 
deviation (RMSD) of 0.3 A (612 Ca pairs). The jRdRp 
shares a high sequence identity with DENY RdRp 
(66%) and WNY RdRp (84%) (Supplementary 
Figure S3) (27,28). The structure of RdRp adopts the 
right-handed topology and shows the presence of the 
three domains characteristic of replicative DNA and 
RNA polymerases, namely, the palm (368-387, 498-545, 
602-716), fingers (274-367, 388^97, 546-601) and thumb 
(717-889) (Figure 2A). The palm domain houses the cata- 
lytic residues D668, D669 and D536. A detailed descrip- 
tion of the jRdRp structure is provided in the 
Supplementary Information. 

Motif F and GTP hinding 

The motif F in the jRdRpoTP structure presented here is 
ordered unhke the case of apo- structures of DENY and 
WNY reported before (Figure 2 and Supplementary 
Figure S2). On superimposition of the structure of 
jRdRpcxp with that of available structures of flaviviral 
RdRps, RMSDs of 2.7 A, 3.8 A and 2.1 A were obtained 
with the structures of DENY apo-RdRp (562 Ca pairs; 
2J7U), WNY apo-RdRp (571 Ca pairs; 2HFZ.pdb) and 
DENY apo-RdRp soaked with GTP (562 Ca pairs; 
2J7W), respectively (28). If the motif F substructure and 
the loop stretching from residues 407^23 (sometimes 
referred to as motif G) are removed from jRdRpoiy struc- 
ture, then these RMSD values come down to 1.6 A, 1.8 A 
and 1.5 A (37). Motif F can be divided into three 
subparts— Motif Fl: M457 to K463; F2: K469 to A470; 
and F3: K471 to 1476. This motif is part of the longest 
stretch of amino acids, which is highly conserved among 
Flaviviruses (454-490) (Supplementary Figure S3). 

The structure of the jRdRpoxp complex shows that 
motif F adopts a beta-hairpin substructure that is 
oriented toward the active site center (Figure 3A and 
Supplementary Figure S4A). The jRdRp^xp complex 
also has a similar motif F substructure. The region just 
beyond the motif F in DENY (474-479) and WNY (475- 
480) forms the first turn of a helix (a 10) and this helix is 
oriented in a direction perpendicular to the active site loop 
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Figure 1. Role of GTP in initiation and elongation by jRdRp. The results of initiation assay in the absence of primer with varying GTP (A) and 
ATP (B) concentrations are displayed. For the initiation assay, 6FAM-labeled UTP was used to visualize the reaction products. Lane 1 in (A) shows 
the migration of molecular markers, and the bands corresponding to 5, 12, 15 and lOmer ssRNA labeled with 6-FAM at the 3' end are marked. Lane 
2 corresponds to the reaction mix with the catalytic mutant of jRdRp, D668A. The wedge and the corresponding values show the increasing 
concentration of GTP (A) or ATP (B) in the reaction mix. (C) Graph showing the cumulative intensity of initiation products obtained in the presence 
of increasing concentrations of GTP and ATP. Error bars represent the standard deviation (independent assays, n = 3). (D) Primer extension assay 
with the jRdRp enzyme. The primer is labeled with fluorophore 6FAM at 5' end as indicated by the cartoon. Lane 1 is reaction mix without Mn^^ 
and Lane 2 is reaction mix with the catalytic mutant, D668A RdRp. The maximally extended primer is marked by an arrow. GTP concentration 
varies across the lanes as indicated by the wedge. ATP, GTP and UTP were maintained at identical concentrations of 100|.iM, each. (E) Graph 
showing the fraction of total primer extended for different amounts of GTP and ATP. Error bars represent the standard deviation (independent 
assays, n = 3). 



(Supplementary Figures S4B and S4C). In the case of 
jRdRp, the equivalent stretch (475-480) adopts an 
extended conformation and is part of the beta-hairpin 
substructure adopted by motif F. 

To ascertain whether motif F becomes ordered on 
binding GTP/ATP, we determined the structure of 
jRdRp in its apo-state (3.65 A). jRdRp was crystallized 
in the space group P2i2i2i with the cell constants 
of a=119.5A, b= 174.0 A, c = 182.7A and 
a = p = y = 90° (Table 1). The complex was refined to a 
final Rfree of 27.3% and Rcryst of 23.4%. There are four 



molecules of RdRp in the asymmetric unit, namely, mol- 
ecules A, B, C and D. In molecules A and B, motif F is 
ordered and has a beta-hairpin substructure similar to that 
in jRdRpcTP structure (Supplementary Figures S5A and 
S5B). In molecules A and B, the residues E466 and F467, 
which are at the tip of motif F loop, form interactions with 
thumb residues of molecules B and A, respectively. These 
interactions stabilize motif F in a conformation similar to 
that observed in the jRdRpoxp structure. In the case of 
molecule C, the beta-strands of the motif F are ordered, 
whereas the electron density for the loop (residues 
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0=3.0 

Figure 2. Structure of the jRdRpGTP complex. The finger, thumb and pahn domains are colored yellow, orange and cyan, respectively. (A) The 
alpha-helices and beta-strands are numbered from the N-terminal of RdRp construct. The side chains of the active site residues (D668, D669 and 
D536) along with GTP are shown in stick representation. (B) The simulated anneahng omit map (Fo-Fc) for GTP is displayed at a contour of 
a = 3.0. 



460-475) connecting the two strands is weak 
(Supplementary Figure S5C). For molecule D, the entire 
motif F is disordered (Figure 3B). The average real space 
correlation for motif F substructure (residue 453^80) is 
0.90, 0.89 and 0.78, respectively, for molecules A, B and C 
in apo-RdRp structure. Unlike for molecules A and B, 
there is no polypeptide chain adjacent to the motif F 
loop of molecules C and D, and consequently, the inter- 
actions formed by E466 and F467 are absent. Because 
there is no external stabilization of motif F in the case 
of molecules C and D, the disorder exhibited by the 
motif F in these molecules would be a closer representa- 
tion of their state in solution. 

The structure of JEV NS5 in the apo-state was recently 
reported (38). This structure (4K6M) showed an ordered 
motif F in a different conformation than that observed in 
the jRdRpcTP structure (Figure 3C). In NS5, the tip of 
motif F forms interactions with the Mtase domain inter- 
action stabilizing motif F in the observed conformation 



(38). Additionally, the authors suggest that — as in the 
case of NS5 from DENV — the Mtase domain can adopt 
different relative orientations with respect to the RdRp 
domain, and therefore the observed stabilization of the 
motif F would be transient (39). Overall, the observations 
made here and the data available for JEV apo-NS5, WNV 
and DENV RdRp suggest that motif F substructure 
exhibits conformational heterogeneity in flaviviruses and 
can adopt multiple structures in absence of interaction 
with NTP. 

Motif F in RdRp (apo- or ligand-bound form) from 
HCV, bovine viral diarrhea virus, norovirus, foot-and- 
mouth disease virus, poliovirus and human immunodefi- 
ciency virus -1 reverse transcriptase adopts a similar beta- 
stranded substructure (40) (Supplementary Figure S6). 
However, in the jRdRp structure, the hairpin loop 
within the motif F substructure points downward and ob- 
structs the tunnel formed between the fingers and thumb 
domain through which NTP enters into the active site 
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(Figure 3D). The length and composhion of motif F in 
these enzymes is different from that of flaviviral RdRps 
and it does not occlude the NTP entry tunnel. 
Additionally, in the case of apo- jRdRp (molecules C 
and D) and apo- NS5 from JEV, the NTP entry tunnel 
is not blocked by motif F. 

Interactions of jRdRp with GTP 

We observed clear electron density for GTP, and in the 
final refined__ structure, the average B-factor of the nucleo- 
tide is 67 (Figure 2B and Table 1). To obtain the 



structure of DENV in complex with GTP, Yap and col- 
leagues soaked apo-RdRp crystals with GTP (28). In the 
consequent structure, only the triphosphate moiety of the 
GTP could be modeled. However, in the jRdRpoxp struc- 
ture, the entire GTP molecule could be placed in the 
observed electron density. 

In the jRdRpoxp structure, the triphosphate moiety of 
GTP interacts with the side chains of basic residues from 
motif F (R460, K463, K471 and R474) and from motif E 
(R734 and R742) (Figure 4A and B). In a recent report, 
Iglesias and colleagues have shown for DENV that the 




Figure 3. Motif F in structures of j'RdRp. (A-C) Stereo images of the G loop (green) and the motif F (red) substructure of jRdRpoTP, apo- jRdRp 
(molecule D) and apo- NS5 from JEV (4K6M), respectively. Residues 398-526 are shown. In apo- jRdRp. the motif F and G-loop are disordered 
(B). In apo- NS5 (C) and jRdRpcxp structure (A), the motif F forms a beta-stranded substructure, but the loop connecting the two strands adopts 
distinct structures. (D) Stereo image of surface representation of the jRdRp showing the NTP tunnel in the jRdRpcxp structure. The polymerase 
surface corresponding to palm, fingers and thumb is colored cyan, yellow and orange, respectively. The RNA binding groove is on backside of the 
displayed surface. Motif F (red-ribbon and surface) occludes the NTP tunnel. 

(continued) 
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Figure 3. Continued. 



residues K463 and K462 are responsible for specifically 
selecting viral template for initiation (41). These residues 
are part of motif F and K463 does interact with GTP. The 
residues of the priming loop (S799 side chain hydroxyl 
group, main chain atoms of W800 and S801) also 
interact with the triphosphate and sugar moieties of 
GTP. The conformation of the priming loop in 
jRdRpoTP is nearly identical to that in DENV RdRp 
structures (apo or GTP soaked) as this stretch superim- 
poses with an RMSD of 0.4 A (26 Ca atoms). In the case 
of apo-RdRp from WNV, the RMSD for superirnposition 
of the priming loop with that of jRdRp is 1.6 A (26 Cot 
atoms). Four residues from motif E — C714, R734, N716 
and Y771 — form a network of water-mediated hydrogen 
bonds with the gamma phosphate of GTP. From motif C, 
the catalytic residue D668 interacts with 2'0 of the ribose 
sugar. D668 forms a water-mediated hydrogen bond with 
the N3 of the guanine base in one of the molecules in the 
asymmetric subunit. 

The carboxyl group of D541 (motif A) interacts with Nl 
(2.7 A) and N2 (3.0 A) of the guanine base. Also, the 06 
atom of guanine base makes a H-bond (3.1 A) with the 
hydroxyl group of S604 (motif B). These interactions with 
the N2 and 06 atoms are specific for guanine base and are 
not possible with the adenine base of ATP. Almost all the 
aforementioned residues that interact with GTP are iden- 
tical in different flaviviral RdRps. The residues S799 and 
N716 are unique to jRdRp and replaced by threonine 
and histidine, respectively, in other flaviviruses 
(Supplementary Figure S3). Overall, the GTP molecule 
is stabilized in the jRdRp active site through interactions 
with residues of the A, B, C, E and F motifs and that of 
the priming loop. 

jRdRpATP 

The backbone structure of jRdRp when bound to GTP or 
ATP is nearly identical (Supplementary Figure S7A), and 
the two structures superimpose with an RMSD value of 



0.2 (612 Ca pairs). The electron density map (simulated 
annealing Fo-Fc omit map) for ATP in the jRdRpAxp 
structure was not as clear as the one observed for GTP 
and exhibited only partial density for the base moiety 
(Figure 5A). This observation suggests that the base 
moiety of ATP is not as well-stabilized as in the case of 
GTP and exhibits rotational flexibility. Comparison of the 
conformation of GTP and ATP in the jRdRpAxp and 
jRdRpcTP complexes shows that the triphosphate region 
of ATP and GTP ahgns well, but the sugar and base 
moieties adopt distinct conformations (Figure 5B). The 
residues of jRdRp that interact with the triphosphate 
moiety in the jRdRpAxp and jRdRpcxp structures are 
identical and attain the same conformation. However, 
the residues of jRdRp that are proximal to the sugar 
moiety of ATP (C714, S715 and S801) are distinct from 
that seen in case of jRdRpcTP- Owing to weak electron 
density for the adenine base, it is difficult to be confident 
about the interactions formed between the base atoms and 
residues of jRdRp. 

The interactions formed between residues of motif F 
and the triphosphate moiety of ATP are similar to that 
seen in the case of jRdRpoTP- Consequently, the motif F 
in the jRdRpAxp structure is ordered and exhibits a con- 
formation similar to that seen in the case of jRdRpoxp- 
As in the case of jRdRpAxp. the motif F occludes the 
NTP entry tunnel in the jRdRpAxp structure also 
(Supplementary Figure S7B). 

Diversity in conformation and binding sites of NTPs 

The location and conformation of bound GTP was 
compared with NTP in other available RdRp-NTP 
complexes (Supplementary Figure S8). In the available 
structure of a flavivirus (DENV) RdRp in complex with 
GTP (2J7W.pdb), the triphosphate moiety could be 
modeled, but this moiety refined with high B-factors 
(>100 A"). In most of the available structures of RdRp- 
NTP complexes, the NTP ligand refines to high B-factors 



Nucleic Acids Research, 2014, Vol. 42, No. 4 2767 




Figure 4. jRdRp-GTP interactions. The finger, thumb and pahii 
domains are colored yellow, orange and cyan, respectively. GTP and 
the interacting amino acid residue are shown in stick representation. 
The water molecules responsible for water-mediated interactions 
between jRdRp and GTP are depicted in the form of red spheres. 
(A) The residues that interact with the triphosphate and sugar 
moieties of GTP are displayed. (B) The residues that form interactions 
with the guanine base and ribose sugar of GTP displayed. The base- 
specific hydrogen bonds formed between D541 and S604 with the N2 
and 06 atoms of guanine are highlighted by dotted lines. 



and/or <1 occupancy. The triphosphate moiety shows 
some overlap with that of GTP in the jRdRpGTP structure 
(Supplementary Figure S8A) (28). In the RdRp-GTP 
structure (lS49.pdb) of bovine viral diarrhea virus 
(family: Flaviviridae; genus: Pestivirus), the base and 
sugar of GTP points to a different direction on compari- 
son (Supplementary Figure S8B) (42). In the case of HCV 
(family: Flaviviridae; genus: Hepacivirus), multiple struc- 
tures of RdRp-GTP or RdRp-UTP complex are available 
(43^5). In these structures, GTP and UTP bind at shghtly 




Figure 5. Conformation of ATP in JRdRpAxp- (A) The simulated 
annealing omit map for ATP is displayed at a contour of a = 3.0. 
(B) The bound conformations of ATP and GTP are compared by 
superimposing the jRdRpo-pp and jRdRpAxp structures (RMSD = 
0.2 A, 612 C-a pairs). 



different locations in distinct conformations that show 
no overlap with GTP in jRdRpoTP (Supplementary 
Figure S8C-S8G). 

The structures of RdRp in coinplex with ATP 
(IHll.pdb) and GTP (2JLG.pdb) for the phi6 bacterio- 
phage (family: Cystoviridae; genus: Cystovirus) are also 
available (46,47). For 2JLG, GTP attains different con- 
formations in the three RdRp molecules in the asymmetric 
unit. In one of the RdRp molecules in the asymmetric 
unit, two molecules of GTP are bound in a location and 
conformation similar to that observed in the full initiation 
complex for this enzyme (lUVN). In the other two mol- 
ecules, the orientation of GTP (and that of ATP in 
IHll.pdb) is anti-parallel to the GTP conformation in 
jRdRpoTP (Supplementary Figure S4H). For Rotavirus 
(family: Reoviridae; genus: Rotavirus) RdRp the GTP 
molecule (2R7X.pdb) attains an inverted orientation at 
a location slightly displaced compared with GTP in 
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jRdRpGTP (Supplementary Figure S4I). The GTP 
molecule is present in a location that overlaps with the 
priming loop in jRdRpoxp. In the RdRp-NTP structures 
from PoHovirus (family: Picornaviridae), foot and mouth 
disease virus (family: Picornaviridae) and norovirus 
(family: Caliciviridae), the conformation of the bound 
NTP is similar to that of the incoming nucleotide during 
the elongation phase (48-50). The binding site and orien- 
tation of the bound NTPs in these structures does not 
overlap with that of GTP in jRdRpoxp (Supplementary 
Figure S8J-S8L). 

Overall, the GTP binding-site location and conform- 
ation in jRdRpoTP was different from that for NTP in 
other known structures of viral RdRp-NTP complexes. 
Thus, GTP binds in the active site of jRdRp in a novel 
orientation. 



Functional relevance of jRdRpcTP structure 

The structures of the initiation state of RdRp molecule 
from phi6 bacteriophage (lUVN.pdb) and reovirus 
(INlH.pdb) are available (46,51). These structures 
contain the RdRp molecule bound to an RNA template, 
the first two 3' nucleotides of which form base pairs with 
two NTPs of the nascent primer in the presence of two 
catalytic ions. The jRdRp structure aligns better with the 
phi6 RdRp structure than the reovirus RdRp enzyme. The 
jRdRpcTP and phi6 RdRpRNA:GTP:GTP structures were 
ahgned through superiniposition of the palm domains. 
As Figure 6 shows, although the palm domain and the 
active site motif aligns well, the GTP molecule in the 
jRdRpoxp structure are oriented almost perpendicular to 
the first two nucleotides (GTPs) present base paired to 



Motif Fp 





Figure 6. Comparison of the structures of jRdRptn-p complex with that of phi6 RcIRprnaiGtpigtp (A) The superimposition of the two structures is 
displayed and tlie C-alplia trace of jRdRp and phi6 RdRp (lUVN.pdb) is colored magenta and green, respectively. The RNA template (colored 
according to element with carbon atoms in yellow) and two GTP molecules (green) from phi6 RdRpRNA:GTP:GTP are displayed in stick represen- 
tation. GTP in the case of jRdRpGxp is colored slate blue. (B) Superposition of the nucleotide and the RNA template is shown separately. The 
comparison shows that the GTP molecule in jRdRp is present in a unique conformation in a distinct binding site. In the full initiation complex of 
jRdRp, GTPj is expected to reorient and adopt the same conformation as GTPp2 for productive catalysis to occur. Also, ATP bound to jRdRp may 
reorient to adopt the conformation exhibited by GTPpi. 
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RNA in the phi6 RcIRprnaigtpigtp complex. We 
identified four residues— D541, R734, R742 and S799— 
that interact with the GTP in jRdRp, but the equivalent 
residues (D329, N317, N324 and N626, respectively) do 
not exhibit functional interaction with template or nucleo- 
tides in the phi6 RdRpRNA:GTP:GTP complex. The residues 
R734 and R742 are part of a substructure spanning from 
residues 733-745 in jRdRpoxp. The ahgnment of 
jRdRpoTP with phi6 RdRpRNA:GTP:GTP shows that this 
region will have to move away from the active site to ac- 
commodate the first incoming nucleotide of the nascent 
primer. 

To ascertain whether the jRdRpcxp structure is 
functionally relevant, the residues D541, R734, R742 
and S799 were mutated to alanine. The residue S799 was 
also mutated to tyrosine. We also generated single and 
double alanine mutants of S604 and D541 — the two 



residues that interact with the Watson-Crick edge of the 
guanine base of GTP. All the mutants could be purified in 
soluble form and were assessed for their abihty to initiate 
RNA synthesis in comparison with the wt enzyme. If the 
observed conformation of GTP in the jRdRpoTP structure 
represents a critical step during initiation, these mutations 
should perturb de novo synthesis without affecting elong- 
ation. We observed almost complete abolition of initiation 
for the D541A, S604A, R734A and R742A mutants 
(Figure 7A and B). 

In the case of S799A and S799Y mutants, the abihty to 
initiate RNA synthesis was enhanced with respect to the 
wt enzyme (Figure 7A and B). The side chain hydroxyl of 
S799 residue makes H-bond with the main chain nitrogen 
of S801 and with the side chain hydroxyl of S804, both of 
which will be disrupted if S799 is mutated to alanine 
(Supplementary Figure S9). S799, S801 and S804 
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Figure 7. Initiation and elongation activity of mutant and wt jRdRp. The results of primer-free initiation assay (A, B) and primer extension assay 
(C, D) for different mutants and wt jRdRp are displayed. The level of initiation and extension activity of different mutants of jRdRp relative to wt 
enzyme are shown in (B) and (D), respectively. Error bars represent the standard deviation (independent assays, n = 3). For the initiation assays (A), 
6FAM-labeled UTP was used. Lane 1 in (A) shows the migration of molecular markers and the bands corresponding to 5, 12, 15 and 20mer ssRNA 
labeled with 6-FAM at the 3'-end are marked. For extension assays, a 20mer primer strand with a 6FAM label at the 5'-end was annealed to a 
complementary 40mer template strand, as indicated by the cartoon in (C). To monitor the extension reaction, a reaction mix without Mn^^ was 
loaded in Lane 1 of the gel displayed in (C) to visualize the unextended primer. Lane 2 in both (A) and (C) correspond to the reaction mix with the 
catalytic mutant of jRdRp, D668A. For the extension assay, the bands corresponding to maximally extended primer (40 nt) in (C) are marked by an 
arrow. 
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residues are part of the priming loop, and this loop has to 
move out of the active site for the double-stranded 
product RNA to extrude out from the enzyme. It is 
possible that the priming loop structure is more mobile 
in the absence of interactions contributed by S799 or 
due to steric interactions formed by the aromatic ring of 
tyrosine. As a result, it might be easier for the product 
RNA to extrude, leading to the observed higher activity. 

Primer extension assays showed that the abihty of wt 
jRdRp, R734A, R742A, S799A and S799Y to catalyze 
elongation was almost identical (Figure 7C and D). In 
the case of the D541A and S604A mutants, the elongation 
activity was marginally reduced. For the double mutant 
S604A-D541A, there was no detectable primer extension 
activity. From these assays, it is evident that the mutations 
have a much more pronounced effect on initiation than on 
elongation. Therefore, the observed bound conformation 
of GTP is mechanistically relevant. 

GTP binding reduces the affinity of jRdRp for RNA in 
the absence of Mn^* 

We assessed binding of jRdRp to a fluorescently labeled 
(5' 6FAM) ssRNA. The strength of this association was 
estimated by measuring change in fluorescence anisotropy 
values on titrating jRdRp against a fixed amount (5nM) 
of lOmer ssRNA. The apparent dissociation constant 
(KdApp) was calculated from the binding curves 
(Figure 8A) assuming a single binding site (see 
Methods). The KjApp value for jRdRp-ssRNA inter- 
action is 39.3 ± 2.3 nM (Hill slope, n = 1.46). Pre-incuba- 
tion of jRdRp at different concentrations with a 10-fold 
molar excess of GTP or ATP led to a substantial drop in 
affinity of jRdRp for RNA (Figure 8A, Supplementary 
Figure SIO and Table 2). In the presence of GTP and 
ATP, there was marginal change in anisotropy with 
increasing concentrations of RdRp, and the data cannot 
be used to fit a binding isotherm and derive binding con- 
stants (Figure 8A). Also, for a fixed concentration of 
jRdRp and RNA, increasing amounts of GTP led to 
gradual decrease in the amount of anisotropy, suggesting 
that the presence of GTP promoted dissociation of RNA 
and jRdRp (Figure 8B). As a control, GTP was titrated in 
the presence of labeled RNA and in the absence of RdRp, 
and no significant change in fluorescence intensity was 
observed (Supplementary Figure Sll). 

The binding affinity of jRdRp toward RNA in the 
presence of GTP/ATP was restored when Mn^^ ion was 
also included in the reaction buffer (Figure 8A). In the 
presence of Mn^^, the KdApp values for jRdRp-ssRNA- 
GTP-Mn^^ and jRdRp-ssRNA-ATP-Mn^^ are 33.9 nM 
(Hifl slope n =1.12) and 40.4 nM (Hill slope n = 1.30), 
respectively (Table 2). These observations show that the 
presence of GTP and ATP inhibits binding of jRdRp to 
RNA, and the presence of divalent manganese ions 
abolishes this antagonistic effect. It was seen that this 
inhibitory effect of GTP/ATP on RNA binding was lost 
in the case of the R734A mutant (Supplementary 
Figure SI 2). This observation is consistent with the infer- 
ence that binding of GTP/ATP in the active site of jRdRp 
inhibits RNA binding in the absence of the catalytic ion. 
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Figure 8. Affinity of j'RdRp for RNA in tlie presence and absence of 
nucleotide. (A) Binding curves for interaction of labeled (6FAM at the 
5'-end) ssRNA with RdRp in the absence and presence of different 
combinations of NTP and Mn^^. (B) The graph displays the effect of 
increasing concentrations of GTP (in the absence of Mn^^) on the as- 
sociation of jRdRp and labeled ssRNA. For both (A) and (B), the error 
bars represent the standard deviation (independent experiment, n = 3). 



Table 2. Dissociation constants for the interaction of jRdRp with 
RNA in the presence/absence of NTPs and Mn^^ 

ssRNA 



Addition of 
GTP along 
with Mn^^ 



Addition of 
ATP along 
with Mn^^ 



KdApp (nM) 
Hill slope, n 
value 



39.3 ± 2.3 
1.46 ± 0.07 
0.9857 



33.9 ± 0.6 
1.12 ± 0.02 
0.9954 



40.4 ± 2.1 
1.30 ± 0.05 
0.9962 



The goodness of fit for the fitted curves displayed in Supplementary 
Figure SIO is shown as values. 



DISCUSSION 

The RdRp molecule from the phi6 bacteriophage is one of 
the best characterized with regard to de novo RNA syn- 
thesis. This enzyme also exhibits a requirement of high 
GTP concentration for optimal de novo synthesis. The 
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first two template nucleotides in three different segments 
(M, S and L) of the negative-stranded phi6 genome are 
CC, CC and CU, respectively. Also, the first two template 
nucleotides in the complementary positive strand for all 
three segments will be UC. Hence, the first or second or 
both incoming nucleotides will be GTP during any initi- 
ation event for replication of the phi6 genome. 

In the case of phi6 RdRp, it has been shown that Mg^^ 
is the divalent ion required for catalysis. According to the 
mechanism presented by Wright and colleagues, the full 
catalytically competent initiation state can be formed in 
diverse ways either by (i) binding of template followed by 
binding of the two NTPs (GTPs) with Mg^^ or (ii) NT? 
binding in a site that is not compatible with catalysis 
followed by reorientation to base pair with template 
RNA and the subsequent entry of another NTP with 
Mg^^ ions or (iii) both NTPs (at sites that are compatible 
with catalysis) with Mg^^ ions (catalytic ion) followed by 
entry of template RNA (52). 

In the vast majority of known flaviviruses, the terminal 
dinucleotide sequence CU is conserved at the 3'-ends of 
both positive and negative strands of the dsRF. Therefore, 
the requirement for high concentrations of GTP is 
surprising because the first incoming nucleotide would 
be ATP. Assuming a similar final initiation state 
complex in the case of flaviviruses as seen for phi6, the 
data presented in this study suggest the following sequence 
of events during initiation in flaviviruses: (i) Apo-RdRp 
wherein motif F can adopt a multitude of conformations 
(Supplementary Figure S13A). (ii) RdRp binds GTP to 
form pre-initiation complex with concomitant ordering 
of motif F. This motif now blocks the NTP tunnel. In 
the observed conformation, GTP binds in a site that 
overlaps with the predicted RNA binding region. Hence, 
GTP binding will prevent entry of RNA into the active 
site (Supplementary Figure S13B). (iii) RdRp binds RNA 
in the presence of Mn^^; concomitantly, the GTP nucleo- 
tide reorients to base pair with C, the second nucleotide in 
the template. In the case of polymerases, the triphosphate 
moiety of the incoming nucleotide is closely involved in 
the octahedral coordination of the catalytic ion along with 
the catalytic residues (53). Hence, it is possible that the 
presence of Mn^^ may trigger reorientation of GTP so 
that its triphosphate moiety can coordinate the catalytic 
ion along with the catalytic residues D668 and D536. The 
reoriented GTP probably attains a conformation wherein 
it will not interfere with RNA binding and the template 
RNA can enter the active site of jRdRp. The GTP 
molecule can now base pair with the second 3' nucleotide 
C of the template. Therefore, addition of Mn^^ could 
abohsh the inhibitory effect of GTP toward RNA 
binding to jRdRp. Because the residues of the motif F 
interact with the triphosphate moiety of GTP, the reorien- 
tation of GTP in an anti-clockwise direction will lever the 
outward movement of motif F. It is also possible that this 
reorientation will lead to a dissociation of GTP and motif 
F. In either case, the motif F will no longer occlude the 
NTP tunnel (Supplementary Figure S13C). (iv) ATP 
diffuses through the now open NTP tunnel and base 
pairs with the first nucleotide U to form the full initiation 
complex (Supplementary Figure SI 3D). Because all the 



components are in place now, productive catalysis to 
form the first phosphodiester bond can occur, leading to 
the synthesis of a dinucleotide primer that can be 
extended. 

In the case of the jRdRp atp complex, ATP will reorient 
to pair with the first nucleotide on template RNA 
(Supplementary Figures S14A and S14B). Base pairing 
of ATP with the first nucleotide U wifl result in signifi- 
cantly less reorientation — especially for the triphosphate 
moiety — as compared with that in the case of GTP 
base pairing with the second template nucleotide 
(Supplementary Figure SMC). Therefore, it is possible 
that motif F conformation wifl not undergo a drastic 
change, and the NTP entry tunnel will still be occluded. 
Consequently, GTP cannot diffuse in and pair with the 
second nucleotide of the template. Hence, the jRdRpAxp 
represents a dead-end complex that cannot lead to pro- 
ductive initiation, and this inference is supported by bio- 
chemical data (Figure 1). Overall, the second incoming 
nucleotide GTP has to bind first to ensure that on 
pairing with the second nucleotide in template RNA in 
the presence of Mn^^ ions, the motif F moves to open 
the NTP entry tunnel and allow diffusion of the first 
incoming nucleotide into the active site to form a full ini- 
tiation complex that wiU lead to productive catalysis. 

The proposed mechanism would prevent stabilization of 
two nucleotides in the RdRp active site in the absence of 
template RNA and thus prevent non-templated RNA 
synthesis. GTP binding ensures that the two nucleotides 
bind temporally before and after RNA binding. The 
proposed chain of events would avoid adventitious and 
wasteful RNA synthesis that could be erroneous, and 
thus prevent the appearance of mutations in the conserved 
terminal sequences. 

A recent study on NS5 from DENV suggests that this 
enzyme has the ability to synthesize AG dinucleotide in 
the absence of template and primer, and the authors 
propose that this ability enforces conservation of the 
terminal sequences (54). However, the amount of di- 
nucleotide product obtained in the presence of template 
was substantially higher. Hence, the mechanism of initi- 
ation proposed above could be the preferred pathway for 
de novo template RNA synthesis by flaviviral RdRp 
enzymes. Also, it is possible that RdRp enzymes from 
different flaviviruses use distinct pathways for initiation. 

The results presented here for JEV RdRp raise the pos- 
sibility that chemical inhibitors that perturb GTP binding 
can compromise initiation of RNA synthesis and thus 
reduce viral titer. OveraU, the study presented here 
provides a basis for the specific recognition of the initiator 
nucleotide GTP by a flaviviral RdRp enzyme and allows 
formulation of a possible mechanism of initiation that 
prevents non-templated RNA synthesis. The observed 
pre-initiation state and the proposed mechanism of initi- 
ation could be the target of therapeutic intervention to 
combat flaviviral infections. 

ACCESSION NUMBERS 

4HDG 4HDH 4MTP. 



2772 Nucleic Acids Research, 2014, Vol. 42, No. 4 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 



ACKNOWLEDGEMENTS 

The authors thank Taruni Roy and Manjula Devi 
Rangappa for help with cloning, expression and purifica- 
tion of the RdRp enzyme. D.T.N, thanks the X-ray dif- 
fraction facility located in the Molecular Biophysics Unit 
of the Indian Institute of Science (funded by Department 
of Biotechnology [DBT] and Department of Science and 
Technology [DST], Government of India) for facilitating 
screening and data collection. D.T.N, acknowledges the 
help rendered by Dr Hassan Belrhah (ESRF) and Dr 
Babu Manjashetty during data collection at the BM-14 
beamline of ESRF. 



FUNDING 

The Department of Atomic Energy (Government of 
India); Ramanujan feUowship from the Department of 
Science and Technology (Government of India) (to 
D.T.N. ); Data collection at the BM14 beamline of 
ESRF (Grenoble, France) was funded by the BM14 
project- a collaboration between DBT (Government of 
India), EMBL and ESRF. Funding for open access 
charge: Ramanujan feUowship from the Department of 
Science and Technology. 

Conflict of interest statement. None declared. 



REFERENCES 

1. Liiidenbacli,B.D., Thiel.H.J. and Rice,C.M. (2007) Flaviviridae: 
The Viruses and Tlieir Replication, 5th edn. Lippincott Williams & 
Wilkins, Philadelphia. 

2. Ghosh,D. and Basu,A. (2009) Japanese encephalitis-a pathological 
and clinical perspective. PLoS Negl. Trap. Dis., 3, e437. 

3. Chambers,!.;., Hahn,C.S., Galler,R. and Rice,C.M. (1990) 
Flavivirus genome organization, expression, and replication. Anmi. 
Rev. Microbiol., 44, 649-688. 

4. Davidson,A.D. (2009) New insights into flavivirus nonstructural 
protein 5. Adv. Virus Res., 74, 41-101. 

5. Egloff.M.P., Benarrocli,D., Selisko,B., Romette,J.L. and 
Canard,B. (2002) An RNA cap (nucleoside-2'-0-)- 
methyltransferase in the flavivirus RNA polymerase NS5: crystal 
structure and functional characterization. EMBO /., 21, 
2757-2768. 

6. Zhou.Y., Ray,D., Zhao,Y., Dong.H., Ren.S., Li,Z., Guo,Y., 
Bernard,K.A., Shi.P.Y. and Li,H. (2007) Structure and function 
of flavivirus NS5 methyltransferase. /. Virol. . 81, 3891-3903. 

7. Koonin.E.V. (1991) The phylogeny of RNA-dependent RNA 
polymerases of positive-strand RNA viruses. /. Gen. Virol., 72, 
2197-2206. 

8. Khromykh,A.A., Kenney,M.T. and Westaway,E.G. (1998) trans- 
Complementation of flavivirus RNA polymerase gene NS5 by 
using kunjin virus replicon-expressing BHK cells. /. Virol., 72, 
7270-7279. 

9. Kapoor.M., Zhang,L.W., Ramachandra,M., Kusukawa,J., 
Ebner,K.E. and Padmanabhan,R. (1995) Association between Ns3 
and Ns5 proteins of dengue virus type-2 in the putative RNA 
replicase is linked to differential phosphorylation of Ns5. J. Biol. 
Oiem., 270, 19100-19106. 

10. UchilP.D., Kumar,A.V. and Satcliidanandam.V. (2006) Nuclear 
localization of flavivirus RNA synthesis in infected cells. /. Virol., 
80, 5451-5464. 



1 1 . Uchil,P.D. and Satchidanandam.V. (2003) Architecture of the 
flaviviral replication complex. Protease, nuclease, and detergents 
reveal encasement within double-layered membrane compartments. 
J. Biol. Chem., 27», 24388-24398. 

12. Welsch.S., Miller,S., Romero-Brey,!., Merz,A., Bleck,C.K., 
Waltlier,P., Fuller,S.D., Antony,C., Krijnse-Locker,J. and 
Bartenschlager,R. (2009) Composition and three-dimensional 
architecture of the dengue virus replication and assembly sites. 
Cell Host Microbe, 5, 365-375. 

13. Chu,P.W. and Westaway.E.G. (1985) Replication strategy of 
kunjin virus: evidence for recycling role of replicative form RNA 
as template in semiconservative and asymmetric replication. 
Virology. 140, 68-79. 

14. Mazzon.M., Jones.M., Davidson,A., Chain, B. and Jacobs, M. 
(2009) Dengue virus NS5 inhibits interferon-alpha signaling by 
blocking signal transducer and activator of transcription 2 
phosphorylation. /. Infect. Dis., 200, 1261-1270. 

15. van Dijk,A.A., Makeyev,E.V. and Bamford,D.H. (2004) Initiation 
of viral RNA-dependent RNA polymerization. /. Gen. Virol., 85, 
1077-1093. 

16. Paul,A.V., van Boom,J.H., Filippov,D. and Wimmer,E. (1998) 
Protein-primed RNA synthesis by purified poliovirus RNA 
polymerase. Nature, 393, 280-284. 

17. Hagen,M., Tiley,L., Chung,T.D. and Krystal,M. (1995) The role 
of template-primer interactions in cleavage and initiation by the 
influenza virus polymerase. /. Gen. Virol., 76, 603-611. 

18. Ackermann,M. and Padmanabhan,R. (2001) De novo synthesis of 
RNA by the dengue virus RNA-dependent RNA polymerase 
exhibits temperature dependence at the initiation but not 
elongation phase. /. Biol. Chem., 276, 39926-39937. 

19. Nomaguchi.M., Ackermann,M., Yon,C., You,S. and 
Padmanabhan.R. (2003) De novo synthesis of negative-strand 
RNA by dengue virus RNA-dependent RNA polymerase in vitro: 
nucleotide, primer, and template parameters. /. Virol., 77, 
8831-8842. 

20. Sehsko,B., Dutartre,H., Guillemot,!. C, Debarnot,C., 
Benarroch,D., Khromykh,A., Despres,P., Egloff,M.P. and 
Canard,B. (2006) Comparative mechanistic studies of de novo 
RNA synthesis by flavivirus RNA-dependent RNA polymerases. 
Virology, 351, 145-158. 

21. Ranjith-Kumar,C.T., GutshaU,L., Kim,M.J., Sarisky,R.T. and 
Kao,C.C. (2002) Requirements for de novo initiation of RNA 
synthesis by recombinant flaviviral RNA-dependent RNA 
polymerases. /. Virol., 76, 12526-12536. 

22. Yi,G.H., Zhang,C.Y., Cao,S., Wu,H.X. and Wang,Y. (2003) 
De novo RNA synthesis by a recombinant classical swine fever 
virus RNA-dependent RNA polymerase. Eur. J. Biochem., 270, 
4952^961. 

23. Makeyev,E.V. and Bamford.D.H. (2000) Replicase activity of 
purified recombinant protein P2 of double-stranded RNA 
bacteriophage phi6. EMBO J., 19, 124-133. 

24. Ballal,R., Cheema,A., Ahmad,W., Rosen,E.M. and Saha,T. (2009) 
Fluorescent oligonucleotides can serve as suitable alternatives to 
radiolabeled oligonucleotides. /. Bioniol. Tech., 20, 190-194. 

25. Otwinowski,Z. and Minor.W. (1997) Processing of X-ray 
diffraction data collected in oscillation mode. Method EnzvmoL, 
276, 307-326. 

26. Potterton,E., Briggs,P., Turkenburg,M. and Dodson,E. (2003) A 
graphical user interface to the CCP4 program suite. Acta. 
Crystallogr. D, 59, 1131-1137. 

27. Malet,H., Egloff.M.P., Sehsko,B., Butcher,R.E., Wright,P.J., 
Roberts,M., Gruez,A., Sulzenbacher,G., Vonrhein,C., Bricogne,G. 
et al. (2007) Crystal structure of the RNA polymerase domain of 
the west nile virus non-structural protein 5. /. Biol. Chem., 282, 
10678-10689. 

28. Yap,T.L., Xu,T., Chen,Y.L., Malet,H., Egloff,M.P., Canard,B., 
Vasudevan,S.G. and Lescar,J. (2007) Crystal structure of the 
dengue virus RNA-dependent RNA polymerase catalytic domain 
at 1.85-angstrom resolution. /. Virol., 81, 4753-4765. 

29. Brunger,A.T., Adams,P.D., Clore,G.M., DeLano,W.L., Gros,P., 
Grosse-Kunstleve,R.W., Jiang,J.S., Kuszewski,J., Nilges,M., 
Pannu,N.S. et al. (1998) Crystallography & NMR system: a new 
software suite for macromolecular structure determination. Acta. 
Crystallogr. D Biol. Crystallogr., 54, 905-921. 



Nucleic Acids Research, 2014, Vol. 42, No. 4 2773 



30. Wmn,M.D., Murshudov.G.N. and Papiz,M.Z. (2003) 
Macromolecular TLS refinement in REFMAC at moderate 
resolutions. Methods Enzymol., 374, 300-321. 

31. Adams,P.D., Afonine,P.V., Bunkoczi,G., Chen.V.B., DavisJ.W., 
Echols,N., HeaddJ.J., Hung,L.W., Kapral,G.J., Grosse- 
Kunstleve.R.W. el al. (2010) PHENIX: a comprehensive Python- 
based system for macromolecular structure solution. Acta. 
Crystallogr. D Biol. Ciy.^tallogr., 66. 213-221. 

32. Roehrl,M.H., WaiigJ.Y. and Wagiier.G. (2004) A general 
framework for development and data analysis of competitive 
high-throughput screens for small-molecule inhibitors of protein- 
protein interactions by fluorescence polarization. Biochemi.stry, 43, 
16056-16066. 

33. Lovett,S.T. (2004) Encoded errors: mutations and rearrangements 
mediated by misalignment at repetitive DNA sequences. Mol. 
Microbiol.. 52, 1243-1253. 

34. Kim,Y.G., Yoo,J.S., Kini.J.H., Kim,C.M. and Oh,J.W. (2007) 
Biochemical characterization of a recombinant Japanese 
encephalitis virus RNA-dependent RNA polymerase. BMC Mol. 
Biol, 8, 59. 

35. Yu.F.. Hasebe,F., Inoue,S., Mathenge,E.G. and Morita,K. (2007) 
Identification and characterization of RNA-dependent RNA 
polymerase activity in recombinant Japanese encephalitis virus 
NS5 protein. Arch. Virol, 152, 1859-1869. 

36. Uchil.P.D. and Satchidanandam,V. (2003) Characterization of 
RNA synthesis, replication mechanism, and in vitro RNA- 
dependent RNA polymerase activity of Japanese encephalitis 
virus. Virology, 307, 358-371. 

37. Gorbalenya,A.E., Pringle.F.M., Zeddam,J.L., Luke,B.T., 
Cameron,C.E., Kalmakoff.J., Hanzlik.T.N., Gordon,K.H. and 
Ward,V.K. (2002) The palm subdoinain-based active site is 
internally permuted in viral RNA-dependent RNA polymerases of 
an ancient lineage. /. Mol. Biol., 324, 47-62. 

38. Lu,G. and Gong.P. (2013) Crystal Structure of the full-length 
Japanese encephalitis virus NS5 reveals a conserved 
methyltransferase-polymerase interface. PLoS Pathog., 9, 
elO03549. 

39. Bussetta,C. and Choi,K.H. (2012) Dengue virus nonstructural 
protein 5 adopts multiple conformations in solution. Biochemistrv, 
51, 5921-5931. 

40. Bruenn,J.A. (2003) A structural and primary sequence comparison 
of the viral RNA-dependent RNA polymerases. Nucleic Acids 
Res., 31, 1821-1829. 

41. Iglesias.N.G., Filomatori,C.V. and Gamarnik,A.V. (2011) The Fl 
motif of dengue virus polymerase NS5 is involved in promoter- 
dependent RNA synthesis. /. Virol., 85, 5745-5756. 

42. Choi,K.H., Groarke.J.M., Young,D.C., Kuhn,R.J., Smith,J.L., 
Pevear,D.C. and Rossmann,M.G. (2004) The structure of the 



RNA-dependent RNA polymerase from bovine viral diarrhea 
virus establishes the role of GTP in de novo initiation. Proc. Natl 
Acad. Sci. USA, 101, 4425^430. 

43. OTarrell,D., Trowbridge,R., Rowlands,D. and Jager,J. (2003) 
Substrate complexes of hepatitis C virus RNA polymerase (HC- 
J4): structural evidence for nucleotide import and de-novo 
initiation. /. Mol. Biol., 326, 1025-1035. 

44. Bressanelh,S., Tomei,L., Rey,F.A. and De Francesco,R. (2002) 
Structural analysis of the hepatitis C virus RNA polymerase in 
complex with ribonucleotides. /. Virol., 76, 3482-3492. 

45. Harrus,D., Ahmed-El-Sayed,N., Simister,P.C., Miller,S., 
Tricoimet,M., Hagedorn,C.JJ., Mahias,K., Rey,F.A., Astier-Gin,T. 
and Bressaneni,S. (2010) Further insights into the roles of GTP 
and the C terminus of the hepatitis C virus polymerase in the 
initiation of RNA synthesis. /. Biol. Chem., 285, 32906-32918. 

46. Butcher,S.J., Grimes,J.M., Makeyev,E.V., Bamford,D.H. and 
Stuart,D.I. (2001) A mechanism for initiating RNA-dependent 
RNA polymerization. Nature, 410, 235-240. 

47. Poranen,M.M., Salgado,P.S., Koivunen,M.R., Wright, S., 
Bamford,D.H., Stuart,D.I. and Grimes,J.M. (2008) Structural 
explanation for the role of Mn^^ in the activity of phi6 RNA- 
dependent RNA polymerase. Nucleic Acids Res., 36, 6633-6644. 

48. Tliompson,A.A. and Peersen,O.B. (2004) Structural basis for 
proteolysis-dependent activation of the poliovirus RNA-dependent 
RNA polymerase. EMBO J., 23, 3462-3471. 

49. Ferrer-Orta,C., Arias,A., Perez-Luque,R., Escarmis,C., 
Domingo, E. and Verdaguer,N. (2007) Sequential structures 
provide insights into the fideUty of RNA replication. Proc. Natl 
Acad. Sci. USA, 104, 9463-9468. 

50. Zamyatkin,D.F., Parra.F., Macliin,A., Grochulski,P. and Ng,K.K. 
(2009) Binding of 2'-amino-2'-deoxycytidine-5'-triphosphate to 
norovirus polymerase induces rearrangement of the active site. /. 
Mol. Biol., 390, 10-16. 

51. Salgado,P.S., Makeyev,E.V., Butcher,S.J., Bamford,D.H., 
Stuart,D.I. and Grimes,J.M. (2004) The structural basis for RNA 
specificity and Ca2+ inhibition of an RNA-dependent RNA 
polymerase. Structure, 12, 307-316. 

52. Wright,S., Poranen,M.M., Bamford,D.H., Stuart,D.I. and 
Grimes.J.M. (2012) Noncatalytic ions direct the RNA-dependent 
RNA polymerase of bacterial double-stranded RNA virus varphi6 
from de novo initiation to elongation. /. Virol., 86, 2837-2849. 

53. Steitz,T.A. (1998) A mechanism for all polymerases. Nature, 391, 
231-232. 

54. SeIisko,B., Potisopon,S., Agred,R., Priet,S., Varlet,!., Thillier,Y., 
Sallamand,C., Debart,F., Vasseur,J.J. and Canard,B. (2012) 
Molecular basis for nucleotide conservation at the ends of the 
dengue virus genome. PLoS Pathog., 8, el002912. 



